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Abstract 

We solve numerically the nonlinear differential equation for the Hele-Shaw, 
Saffman- Taylor problem derived in the preceding work. Stationary solutions 
with no free phenomenological parameters are found to fit the measured pat- 
terns. The calculated finger half-widths as a function of the physical parameters 
of the cell, compare satisfactorily with experiment. PACS numbers: 47.20.Dr, 
47.54.+r, 68.10.-m 



1 Stationary finger solutions in the Hele-Shaw cell 



In the preceding work[l], we developed nonlinear differential equations for the finger- 
ing phenomenon in the Hele-Shaw cell. [2] 

The time independent equations for the stationary finger are 
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where the suffix indicates differentiation with respect to x, rj is the y coordinate of 

the finger profile. — = — — , with ii the viscosity of the displaced fluid, T, the 
B T b 2 

effective surface tension and b the thickness of the Hele-Shaw cell.[l] Equation (1), 
expressed in terms of the arclength measured from the tail of the finger becomes 

(2) 

V 

In eqs.(l,2), e = — . With U the velocity of the tip of the nose of the finger in 

the laboratory (cell stationary) frame, and, V a velocity parameter determined self- 
consistently from the existence of a stationary solution. Asymptotically far back at the 
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tail we must have W(rj) = 0, or, e = — — . The parameter e, is then determined 

by the solutions to the equations and is not a phenomenological parameter. 

Equation (1) can be transformed to a second order one in terms of the angle 
tangent to the curve, with rj obtained by integration. For ds starting at the tail, 
where 6 ~ ty, ds cos(6) = — dx. 

Equation (2) reads 

r\ = A — J ds sin{6) (3) 
With W{rj) defined in eq.(l). 

In the preceding work we analysed qualitatively the predicted finger half-widths 
A, as a function of the parameter B. We found e = 1 in the limit of B — >• oo. In this 
limit the potential determines the finger half-width to be A = 1, as expected. At the 
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other end of B — > 0, or infinite viscosity limit, we argued [1], that e = in the spirit 
of the no-slip condition for viscous fluids. In this limit the potential W(rj), yields 

A = -j= « 0.447. Figure 8 of Tabeling et al.[6], shows that for — — >• oo the finger 

half-width for various aspect ratios tends to A ~ 0.45. 

Therefore we can be quite confident on the potential W. On the other hand, the 
dynamics as described by the third order differential equation linear in the parameter 
B is certainly limited to low capillary numbers for which only the lowest order Young- 
Laplace expression suffices. At larger capillary numbers higher order terms in the 
curvature are needed. [7] In the following we limit our calculations to relatively low 
capillary numbers. 

The numerical solution was started from positive y = A and followed around the 
finger to y = —A. We employed a fourth order Runge-Kutta algorithm. We opted 
for the mixed integrodifferential approach of eq.(3), that proved easier to handle 
numerically. 

For fixed B we varied the value of the asymptotic half- width (thereby fixing e), 
with 9'(s =) = and 9(s = 0) = 180°, until a solution is found whose nose tip shows 
up at exactly 90°. Moreover, only solutions where the angle decreases smoothly from 
the tail forwards were accepted. We rejected meandering (even slightly so) solutions 
that could appear as a many pronged finger. 

We proceed to show 3 profiles found for finger half-widths A = 0.93, 0.82, 0.66. 
They are compared with the analytical solution of Pitts [5] that fits the data extremely 
well for finger half-widths smaller than around A ~ 0.8 and lies slightly inside the 
data points for higher values of A. The figures are depicted for rescaled coordinates 
in terms of A as has became customary [3, 5]. 

From figures 2-4 it is clear that the solutions match the phenomenological solutions 

of Pitts cos( — — )e~2~ = 1, and are even superior to it for large values of A. A 
2 A 

graph of the asymptotic finger half-width A as a function of B, is shown in figure 
5. The present results improve upon the calculations of McLean and Saffman[3] for 
small capillary numbers. Note however, that our results do not include finite gap 
corrections that lift the curve upwards for small B, whereas the results of McLean 
and Saffman[3]. shown in the figure do adjust for such corrections. Without these 
corrections the latter results fall much below the depicted curve. 

We found numerically stationary finger profiles for B < 50. Beyond this value, 
the large contributions introduced by the potential in eq.(3), render the integration 

from the tail unstable. The curve of A versus — Several attempts were made in the 

present study, to improve the algorithm by matching solutions run from the tail and 
from the nose. However, starting from the nose is very problematic, because of the 
singularity in the slope. 

We found that for a a fixed value of B, there corresponds a single value of A; a 
unique finger solution. The selection of the finger half-width is univocous. However, 
we cannot rule out completely the existence of other unstable branches not accessed by 
the numerical method, as found by Vanden-Broek[4] for the McLean and Saffman[3] 
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Figure 1: Reseated finger profile for X = 0.93 at ^ — 12. Numerical solution, solid 
line, analytical solution of Pitts [10], dashed line. 
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Figure 2: Reseated finger profile for X = 0.82 at ^ — 20. Numerical solution, solid 
line, analytical solution of Pitts [10], dashed line. 
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Figure 3: Rescaled finger profile for A = 0.66 at ^ — 38. Numerical solution, solid 
line, analytical solution of Pitts [10], dashed line. 
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Figure 4: Asymptotic half-width of the finger, A as a function of j^. Present work, 
full line, theoretical results of McLean and Saffman [7] corrected for finite film effects, 
dashed line, and, experimental results of Tabeling et al. [9], dash-dot line. 
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set of integrodifferential equations. 



2 Conclusion 

The numerical stationary profiles obtained from the nonlinear differential equation 
derived in [1] , fit well the experimental ones, that are, in turn, quite accurately repro- 
duced by the phenomenological solution of Pitts[5]. The predicted finger half-width 
A, as a function of capillary number, fits the experimental results in the low capillary 
number regime. In the course of the numerical investigation we found solutions that 
correspond to various multiple finger-like structures. These structures originate from 
specific boundary conditions at y — A. The existence of such complicated patterns is 
quite encouraging, especially if we want to proceed further to the production, devel- 
opment, competition and splitting between fingers with the time-dependent equation 
displayed above. 
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